Single exciton spectroscopy of single-Mn doped InAs quantum dots 
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The optical spectroscopy of a single InAs quantum dot doped with a single Mn atom is studied 
using a model Hamiltonian that includes the exchange interactions between the spins of the quan- 
tum dot electron-hole pair, the Mn atom and the acceptor hole. Our model permits to link the 
photoluminescence spectra to the Mn spin states after photon emission. We focus on the relation 
between the charge state of the Mn, A or A~ , and the different spectra which result through either 
band-to-band or band-to-acceptor transitions. We consider both neutral and negatively charged 
dots. Our model is able to account for recent experimental results on single Mn doped InAs PL 
spectra and can be used to account for future experiments in GaAs quantum dots. Similarities and 
differences with the case of single Mn doped CdTe quantum dots are discussed. 

PACS numbers: 



I. INTRODUCTION 

Probing a single magnetic atom in a solid state en- 
vironement is now possible by scanning tunneling mi- 
croscopy (STM), both in metallic^ and semiconducting 
surfaces 2 ^^, and by single exciton spectroscopy in semi- 
conductor quantum dotsSiLS, among other techniques. 
These experiments permit to address a single quantum 
object, the spin of the magnetic atom, and to study its 
exchange interactions with surrounding carriers. Quan- 
tum dots doped with a single magnetic atom are a model 
system for nanospintronic o 9 : 10 ' 1 1 . These systems can also 
shed light on the issue of whether or not it is possible to 
dope semiconductor nanocrystals^. 

The focus of this work is the single exciton spec- 
troscopy of a single Mn atom in a InAs quantum dot 
(QD), motivated by recent experimental results on InAs 
QD^ and keeping in mind the relation to previous ex- 
periments on single Mn doped CdTe£. The photolumi- 
nescence (PL) spectra of a CdTe QD doped with only 
one Mn atom display six narrow peaks, each of which 
corresponds^ to one of the six quantum states of the 
5 = 5/2 multiplet formed by the 5 Mn d electrons^ in 
Mn 2+ . Hence, the spin state of the single Mn atom af- 
ter its interaction with an exciton in a CdTe QD can 
be read from the energy and polarization of the emit- 
ted photonSH. Another interesting result in CdTe dots 
doped with Mn, is the fact that the PL spectrum changes 
radically when a single carrier is electrically injected into 
the doti This has been explained in terms of the dif- 
ferent effective spin Hamiltonian for the Mn as a single 
additional carrier is added into the doli^ ' 11 ! 15 . 

Since Mn is an acceptor in InAs quantum dot, the sin- 
gle Mn exciton PL is expected to be different from CdTe, 
where Mn acts as a isoelectronic impurity. Even before 
photoexcitation, charge neutrality implies that the Mn 
acceptor in InAs binds a hole. In this neutral acceptor 
complex, A , the spin of the Mn is antiferromagnctically 
coupled to the spin of the acceptor hole, so that A be- 



haves as effective spin F = 1 object. When a neutral 
exciton X° is created in a InAs dot doped with 1 Mn, 
there are 4 spins interacting: the QD electron, the QD 
hole, the Mn and the acceptor hole. Indeed, recent ex- 
perimental observations report a band-to-band transition 
(X°A° — > A ) PL spectrum with 5 peaks with different 
intensities at zero applied field, instead of the 6 almost 
identical peaks in CdTe. The presence of the acceptor 
hole also opens a new optical recombination channel, the 
band-to-acceptor transition (X°A° — ► h + A~) such that 
the conduction band electron ionizes the Mn acceptor 
without filling the quantum dot hole. 
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FIG. 1: Schematic energy levels of the InAs dot doped with 1 
Mn. Upper row: band-to-band transition. Lower row: band- 
to-acceptor transition. (Color online) 

The goal of this paper is to provide a theoretical frame- 
work to understand the relation between PL spectra of 
a single InAs quantum dot doped with one Mn and the 
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interactions, charge and spin state of the relevant degrees 
of freedom. The rest of this paper is organized as follows. 
In section II we describe the theoretical framework, in- 
cluding the spin models for the relevant degrees of free- 
dom, and the framework to calculate PL. We discuss both 
a 4-spin-model and a simpler 2-spin-model, proposed in 
reference H, and how they are related. In section III we 
present our simulations for the PL of neutral quantum 
dots. We consider both band-to-band transitions, such 
that the final state is the neutral acceptor A , and the 
band-to-acceptor transition, such that the final state is 
h + A~, i.e. a QD hole interacting with the spin S = 5/2 
of the ionized Mn acceptor. We find that an antifcr- 
romagnetic QD hole-Mn coupling can still yield an ef- 
fective ferromagnetic coupling between the QD hole and 
the Mn-acceptor complex (consisting on the Mn ion and 
the acceptor hole), as observed by Kudclski et aL— . In 
section IV we present our results for negatively charged 
quantum dots. In this case, since there are two electrons 
and two holes, there are 2 recombination pathways and 
4 possible sets of PL spectra are possible. In section V 
we summarize the main results. 



II. THEORETICAL FRAMEWORK 
A. Photoluminescence and eigenstates 

Our goal is to extract information of the quantum state 
of the single Mn spin in the quantum dot from optical 
spectroscopy data. We adopt a phenomenological ap- 
proach where only the spin degrees of freedom of the Mn, 
quantum dot (QD) carriers and acceptor hole are consid- 
ered. The various spin couplings are chosen to respect the 
symmetries of the problem and are fitted to experimen- 
tal data when available. The relevant electronic states of 
the quantum dot are described in terms of few-spin quan- 
tum states which depend on the different spin exchange 
interactions. 

In the calculation of the PL, it is convenient to dis- 
tinguish between the ground state and the exciton state 
manifolds^. Both the number of states and the effec- 
tive spin Hamiltonian of these manifolds depend on the 
charge of the dol£ii and on the charge state of the ac- 
ceptor complex. For instance, the relevant degrees of 
freedom of the ground state manifold (GSM) of a neu- 
tral QD are the spins of the Mn spin and the acceptor 
hole. The exciton state manifold (XSM) enlarges the 
GSM with the addition of the QD electron and the QD 
hole. By definition, the GSM is defined by the eigenstates 
of the Hamiltonian of the dot before the photoexciton is 
injected: 

H G \^ G ) = E G \V G ) (1) 

In the case of Mn in a neutral InAs QD, G runs over 
the 24 possible states that can be formed with a spin 
S = 5/2 of the Mn and a pseudospin 3/2 of the accep- 
tor. Notice that Eq takes different values (ground state 



spin splittings) due to cither exchange coupling between 
the Mn and the acceptor hole or, in the case of compen- 
sated impurities or II- VI semiconductors, coupling to an 
external field. The XSM is formed by the eigenstates 
of the exciton Hamiltonian, which can be written as the 
sum of Tic and the new terms involving all the couplings 
of the photocarriers with the degrees of freedom before 
excitation: 

Hx\*x)=E x \^ x ) (2) 

Both Ex and Eq and their wavefuntions are obtained 
from diagonalization of the model Hamiltonians de- 
scribed below in detail. 

The PL spectrum for a given polarization state A is 
related to these states^: 

h{uj) = nx\(^ G \Vx\^x)\ 2 S [hw - {E x - E G )] (3) 

X,G 

Here nx is the probability that a given XSM state is oc- 
cupied and (ty g\P x) are the the matrix elements, of 
the interband electric dipole operator— that promotes an 
electron from the valence states to the conduction states 
and viceversa. This operator obeys the standard optical 
selection rules associated to the photon with polarization 
A and does not affect the Mn spin state. Explicit expres- 
sions of this operator are provided once we discuss the 
nature of Hg and Tix and their eigenvectors. 

In a non-magnetic dot, the PL spectrum has a single 
line at zero magnetic field. A distinctive feature of mag- 
netically doped dots is the appearance of several lines at 
zero magnetic field£i£. According to equation the 
appearance of several lines in the PL spectra can occur 
both due to splittings in the GSM and in the XSM. The 
intensity of the lines depends on two factors, the quan- 
tum mechanical matrix elements and the statistical occu- 
pation of the emitting state, nx- This quantity depends 
on the complicated non-equilibrium kinetics of the pho- 
toinjected carriers. Instead of solving a non-equilibrium 
master equation^, we assume that the emitting states 
are in a thermal equilibrium with an effective temper- 
ature which can be larger than the temperature of the 
lattice. This phenomenological approach is supported by 
experimental results in the case of single Mn in a CdTe 
QE£. 

B. 4-spin-model 

1. Ground State Manifold 

Mn has two s electrons which participate in the sp- 
bonding, whereas In has three. Therefore, substitu- 
tional Mn in InAs behaves like an acceptor. EPR 16 
and photoemissionii experiments indicate that Mn re- 
tains the five d electrons when doping concentrations are 
small. Hence, Mn keeps an oxidation state of +2 result- 
ing in an effective charge of —1, that repels the electrons 
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nearby. The Mn impurity remains charge neutral, at the 
scale of a few unit cells, by binding a hole. The binding 
energy of the hole is of 110 meV in Ga(Mn)A o 18 ' 19 and 
28 meV in In(Mn)As^i2£. The acceptor hole state has a 
radius of approximately 1 nm and has been probed by 
STM experiments both in GaAs^ and InAs^. In bulk the 
acceptor hole has a fourfold degeneracy inherited from 
the top of the valence band, which is lifted by quantum 
confinement and/or strain. Because of the strong spin- 
orbit interaction, it is convenient to treat the acceptor 
hole as a spin j = 3/2 object, exchange coupled to the 
Mn spin M = 5/2. 

The operators acting upon this object are the four 
by four J = 3/2 angular momentum matrices, j. In 
the spherical approximation^! the Mn-acceptor hole spin 
coupling reads^: 



H M ,j =eM ■ j 



(4) 



where e = +5 meV is the antiferromagnetic coupling be- 
tween the acceptor hole and the Mn, M are the 5 = 5/2 
spin matrices of the Mn and j are the J = 3/2 matrices 
corresponding to the total angular momentum of the va- 
lence band states. The Hamiltonian equation (0J is read- 
ily diagonalizcd in the basis of the total spin F = M + J, 
the spin of the Mn plus acceptor hole complex. F can 
take integer values between 1 and 4. The eigenvalues are 
E(F) = ^F(F + 1) + -Bo- Since the coupling is antifer- 
romagnetic, the ground state has F = 1, separated from 
the F = 2 states by a relatively large energy barrier of 
2e. Hence, as long as the Mn-acceptor hole complex is 
not distorted by perturbations that couple different F 
manifolds and temperature is low enough, it is a good 
approximation to think of it as being a composite object 
with total spin F = 1. Since they play an important 
role, we explicitly write down the 3 wave functions of the 
F = 1 manifold: 
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where we use the notation 



\F = 1,F Z = ±1,0) = Y, C M ^AF,F z )\M z ,j z ) 

We see that because of the very strong exchange interac- 
tion between the acceptor hole and the Mn spins, their 



spins are strongly correlated and they are not good quan- 
tum numbers separately. 

Following Govorov^, we assume that the QD is larger 
than the bulk acceptor state. Thus, the QD perturbs 
weakly the acceptor state. In this approximation we can 
distinguish between quantum confined or QD states and 
acceptor states. The former are extended all over the 
dot the latter arc tightly bound to the Mn impurity and 
their energy lies in the gap (see figure []}. The opposite 
scenario, in which the QD size is comparable or smaller 
than the acceptor state, has been considered by Climente 
et al<21, would yield different results incompatible with 
the experiments of Kudelski et as discussed below. 

The cubic symmetry of the ideal crystal and the pres- 
ence of quantum confinement and strain result in addi- 
tional terms in the Hamiltonian, that need to be summed 
to equation (|4|). Both quantum dot confinement and 
strain can result in a splitting of the light and heavy 
hole bands, that can be modelled a —Dj z term. This 
term would be present in thin film layers with strain and 
still preserves rotational invariance in the xy plane. The 
presence of the quantum dot potential will break this in 
plane symmetry. To lowest order— this can be modelled 
by an additional term in the Hamiltonian, E(j%. — jy). 
Notice that we assume that these perturbations act on 
the acceptor state only and not on the Mn d electrons. 
This is justified since the hole is spread over tens of unit 
cells, whereas the Mn d states are confined within an unit 
cell. Hence, we take the following model for the ground 
state Hamiltonian: 



H G = eM ■ j - Djl + E(j 2 x 



(6) 



As we discuss below we have e ^> D 3> E. Within the 
F = 1 lowest energy manifold the D term splits the 
triplet into a F z = ±1 doublet and a F z = singlet. 
The E term hybridizes the F z = ±1 states, resulting in 
a small hybridization splitting. 

We obtain the the eigenstates of Hg by expressing 
them as linear combinations of \M z ,j z ) = \M Z ) <g> \ j z ): 



I*g) = 



E 



C G 



\M z ,j z ) 



E 

F,F Z 



V% F JF,F Z ) (7) 



and diagonalizing numerically the Hamiltonian matrix. 
The use of the \F,F Z ) basis might be better for interpre- 
tation of the results. 

Finally, in some instances we need to consider the ion- 
ized acceptor complex, h + A~. This is the case if we 
consider the band-to-acceptor transition in neutral dots 
or if we consider a charged quantum dot. The spin of 
the Mn inside the A~ state is M = 5/2 and should have 
properties similar to those of Mn in CdTo® . 



2. Exciton state manifold 

We now consider states with an electron and a hole 
in the QD lowest energy levels in the conduction and 
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valence band respectively (see figure [!} . In contrast to 
the case of neutral Mn in II- VI semiconductor, the ex- 
citon states involves 4 spins instead of three: the Mn 
(M = 5/2), the QD conduction electron a c = ±1/2, the 
QD hole (T/j = ±3/2 and the acceptor hole J = 3/2. Since 
we ignore LH-HH mixing for the QD valence states, the 
QD holes arc heavy hole, with well defined J z = ±3/2 
(or f|",4)- As a result, the spin couplings of the QD to 
the other spins (Mn, acceptor hole and QD electron) are 
Ising like. Including the small LH-HH mixing present in 
the QD hole state results in a small spin-flip terms in 
the exchange Hamiltonian of the QD holeii. We label 
the hole spin states as the time reversed states of the va- 
lence electronic Bloch states with quantum number cr^&, 
ah = —a v . With this notation, the spin of a given state 
that features one quasiparticle in the valence band, ei- 
ther one electron or one hole, is the same as the spin of 
the quasiparticle. With this notation, the exciton spin 
X satisfies the rule X = 0) x + a c and takes values ±1 
states for optically active excitons and ±2 for optically 
dark exciton o 26 i 27 . 

Since we have 4 spin degrees of freedom we need to 
consider the 6 two spin couplings between them: 

7~Lx = 7~Lg + 7~i c ,h + 7~ih,M + 7~ih,j + 7~ic,M + 7~ic,j + Hz(8) 

One of them, the M ■ j term, is present both in the 
GSM and XSM. The symmetry and the coupling strength 
characterize a given spin-spin interaction. In spin ro- 
tational invariant systems two spins s*i and 5*2 interact 
via Hciscnberg coupling, s\ • S2- When the interplay of 
spin-orbit coupling and lack of spherical symmetry break 
spin rotational symmetry, spins are coupled with differ- 
ent strengths along different directions. An extreme case 
are flat self-assembled quantum dots, for which the low- 
est energy hole states are purely heavy holes such that 
in-plane couplings are strictly forbidden^, resulting in 
Ising couplings. In the opposite limit, the conduction 
band states and the Mn d states have no orbital momen- 
tum which greatly reduces the size of spin-orbit inter- 
actions, resulting in Heisenberg couplings between each 
other. An intermediate situation would be that of holes 
in spherical nanocrystals, where, in spite of strong spin 
orbit interactions, Mn-holc exchange is still described 
with a Heisenberg coupling 2 ^. Following previous work in 
CdTeii^, we take the QD hole-Mn and the QD elcctron- 
Mn coupling as antifcrromagnctic Ising and ferromag- 
netic Heisenberg respectively. 

The second term in equation (jHJ) is the longitudinal QD 
electron- hole exchange that splits the bright ±1 and dark 
±2 excitons in two doublets: 

H e h = +Jeh&c&h (9) 

Since the dark doublet, for which a c &h > is lower in 
energy we have J e h < 0. For simplicity, we neglect trans- 
verse electron-hole exchange. 

In the case of A we also need to include the coupling of 
the QD electron and QD hole spins both to the Mn spin 



and to the acceptor hole. We assume the same symmetry 
for the two couplings: 

7~(hM + 7~(hj = JhMOh,M z + Jhj&hjz (10) 

Notice that the sign of the hole-hole coupling is not clear 
a priori. The QD conduction electron couplings are: 

yicM + ticj = —JcM&c ■ M - JcjVc ■ j (11) 

Finally, we include the Zeeman coupling of the various 
spins to an external magnetic field. For simplicity, we 
ignore the orbital coupling to the magnetic field. In the 
neutral dot there are 2 spins in the GSM and 4 spins in 
the XSM that are coupled to the magnetic field. Thus, 
we need 4 (^-matrices. In this paper we only consider 
magnetic fields along the growth axis. The couplings 
read: 

V.Z = HbB z (g c a c + g h a h + g M M z + g j z ) (12) 

where fj, B = = +0.0579 meV/T. 

The exciton states in the 4-spin-modcl are obtained 
by numerical diagonalization of the Hamiltonian. We 
express them as linear combinations of the product basis 
\M z ,j z ,(T c ,a h ): 

l*x)= E C M^,a c ,a h \M z ,j z ,a c ,a h ) (13) 

M z ,j z ,a c ,cr h 

The 4-spin-model has 96 eigenstates, as many as the 
product of the (25 + 1) x (2 J + 1) x 4, where 2S + 1 = 6 
is the multiplicity of the Mn spin, 2 J + 1 = 4 is the mul- 
tiplicity of the acceptor hole spin, and 4 is the number 
possible of quantum dot exciton states. 

The eigenvalues and eigenvectors of Hx depend on the 
strength of the spin couplings. These depend both on the 
material and the sample. For instance, exchange coupling 
between the QD carriers to the Mn spin is given by the 
product of the material exchange integrals^, a and (3, 
and the probability amplitude of the QD envelope func- 
tions for cither electrons or hole s 13 i 15 , which is clearly 
a sample dependent property. For the same reason, the 
hole-Mn coupling is much stronger for the acceptor state 
than for the QD state. Here we choose the numerical 
values of the exchange coupling constants to account for 
the experimental data. Sample to sample variations will 
result in different PL spectra. 

Importantly, as long as the Mn-acceptor hole is the 
dominant coupling, there are 12 lowest energy exciton 
states, well separated from the rest. These 12 states cor- 
respond to the possible combinations of the 4 exciton 
states and the 3 F = 1 states. Although these states are 
predominantly F = 1, they are somewhat mixed with 
higher F states. It must be noted that, assuming a ther- 
mal occupation with an effective temperature, the PL is 
predominantly given by transitions from the F = 1 man- 
ifold. 
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C. 2-spin- model 

The 4-spin-model has as much as 6 exchange constants 
which might not be possible to extract from comparison 
with PL experiments. The situation can be significantly 
simplified trading off some accuracy. In the GSM, we can 
remove the F > 1 states as long as they are not thermally 
occupied (fc^T < 2e) and not mixed dynamically through 
the terms that break rotational symmetry, i.e. as long as 
both D and E are also much smaller than 2e. If these 
conditions are met, we can use a single spin model for 
the ground stated. The 3 by 3 Hamiltonian in the F = 1 
subspace reads: 



H G = -VFl + £{Fl - Fy) + g Ff i B F z B z 



(14) 



A Hamiltonian similar to this has been used to model Mn 
in GaAs quantum wells^. The coupling constants of the 
2-spin-model are obtained from those of the 4-spin-model 
by representing equations j6]) and (fT2|) in the basis set of 
the F = 1 states from equation ([5]). We obtain 



V 
£ 
ill- 



= m D 

7 

= -t9m 



\9j 



(15) 



The V term splits the F = 1 triplet into a singlet |1, 0) 
and a doublet, |1, ±1). The £ term mixes the two states 
in the ±1 doublet resulting in a bonding and anti-bonding 
states along the Y and X axis respectively. 

Along the same lines, the XSM Hamiltonian can be 
approximated by a simpler 2-spin-model one, if we treat 
the optically active exciton as a quantum Ising degree of 
freedom, X z = ±1, coupled to the spin F = 1 formed 
by the Mn spin and the acceptor hole. In that case, the 
XSM Hamiltonian reads: 



H x =n G + JF Z X Z + g xl i B X z B z 



(16) 



The strength of the effective exciton Mn complex cou- 
pling is related to the bare coupling constants of the 4- 
spin-model through: 



21 9 
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(17) 



In obtaining equation (jTTJ) we set the electron-Mn com- 
plex interaction to zero. It is worth noting that the sign 
of <7f, the effective g-factor of the effective composite 
spin F, as well as the sign of J , the effective exciton- 
Mn complex coupling, could be different from those of 
the constituent particles. In particular, even if both the 
QD hole-acceptor hole J^j and the QD hole-Mn couplings 
JhM are antifcrromagnetic, we could have a negative (fer- 



reversal occurs because the very strong Mn-acceptor hole 
interaction distorts their wave functions and perturbs 
their couplings to a third spin. 
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romagnctic) effective coupling if J^j > 
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weak lines due to mixing by £ {F% — Fy ) 

FIG. 2: Energy level diagram for the neutral exciton band- 
to-band transitions. Both the ground state and exiton state 
(X = +1) energy lines are shown, as well as their evolution 
as a function of the magnetic field. (Color online) 



The advantage of the 2-spin-model is that it can be 
solved analytically. The details are provided in the ap- 
pendix. Importantly, the F z = state is decoupled from 
the F z = ±l pair, both in the GSM and the XSM. The 
energy level diagram is shown in the figure |J2J). The GSM 
features a weakly split doublet. The energy separation 
with the higher energy F z = singlet is T>. The doublet, 
denoted x and y, is a linear combination of the F z = ±1 
states. The small splitting within the doublet is aproxi- 
mately equal to £. The Zeeman term affects the F z = ±1 
states. 

Since in the 2-spin-model the F z = state is decoupled 
from the F z = ±1 states, it is convenient to think of 
the doublet \F = 1,F Z = ±1) as a isospin 1/2 space. 
Both the exchange coupling and the applied field (in the 
Faraday geometry) act as effective magnetic fields along 
the isospin z axis, whereas the £ term acts an effective 
magnetic field in the isospin xy plane. In the 2-spin- 
model the bright exciton does not shift the F z = state. 
Thus, the exciton exchange and the magnetic field mix 
the x and y wave functions of the F z — ±1 doublet. The 
mixing opens two new optical transitions, marked in the 
diagram of figure ([2|). 



D. Optical selection rules 

We now discuss the selection rules associated to the 
dipole Vx operators that we need to use in equation f3]). 
They promote an electron from the valence band to the 
conduction band and viceversa. In a Mn doped III-V 
QD there are two relevant valence band levels so that 
exciton recombination can occur through two channels, 
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as shown in figurc[T] band-to-band and band-to-acceptor. 
These transitions have different energies and, since the 
wavcfunctions of these holes are not the same, different 
optical selection rules. 

Emission of a a + (er _ ) photon in the direction normal 
to the QD layer takes away (adds) one unit of angular 
momentum L z from the system. In the band-to-band 
transitions we shall ignore LH-HH mixing of the quan- 
tum dot hole state. Therefore, emission of a cr + (cr~) 
implies the removal of the +1 (—1) exciton. Hence, in 
the band-to-band channel, the transition operators are 
defined from its action upon XSM states as a projector: 

M z ,j z 

-P+l^x) = £ C^ Lt \M Z7 j z ) (18) 

M z , 3z 

We omit the prefactor proportional to the single parti- 
cle dipole matrix element, which is a convolution of the 
atomic and the envelope wave functions. Ignoring LH- 
HH mixing in the band-to-band transition implies that, 
in our model, linear polarization can only occur through 
quantum coherence between the +1 and —1 exciton. 

In the band-to-acceptor transition the dipole opera- 
tor moves an electron from the QD conduction level to 
the acceptor level, for which we can not ignore LH-HH 
mixing. Hence the gain (loss) of one unit of angular 
momentum upon er + (c _ ) photon emission can occur 
also through the light hole channel. Thus, in a band-to- 
acceptor transition the spin of the annihilated conduction 
band electron a c and the acceptor hole j z are given to the 
a photon: 

a c +j z = ±l (19) 

The band-to-acceptor transition operator V± is fully de- 
scribed by its action upon a given state |\&x) of the XSM. 
The resulting GSM state read, for <j + emission: 

V+\*x) = J2 C M„ +3 /2,i,* h \M z ) ® \a h ) + 

M z ,a h 

+^ ii+1/2iTiff jM i )® \a h ) (20) 
and for o~ emission: 

V-\^x) = C M z ,-3/2,t J M z) ® K> + 

M z ,cr h 

+ 4^,-1/2,1, J Af*> ® Wh) (21) 

Notice that these operators leave the quantum dot hole 
unchanged and connect states in which the Mn spin is 
strongly coupled to the acceptor hole to states where the 
acceptor hole is compensated and the Mn spin is only 
coupled to the QD hole. Notice that, both in band-to- 
band and band-to-acceptor transitions, the dipole oper- 
ators do not act on the Mn d electrons. Therefore, the 



spin of the Mn is conserved during the photon emission 
processes. This is in contrast to intra-shell transitions rel- 
evant when the bandgap is larger than the intra-atomic 
transitions^. 



III. NEUTRAL EXCITON SPECTROSCOPY 

We now present our calculations for the PL spectra 
for neutral InAs quantum dots doped with one Mn. We 
consider both band-to-band X°A° — > A and band-to- 
acccptor X°A° — > h + A~ transitions. The symmetry of 
the state left behind after photon emission is very differ- 
ent in these two cases. In reference HI the band-to-band 
transition X°A° — > A has been experimentally observed 
and described with a model very similar to the one pre- 
sented in the previous section. We first revisit this case 
which permits to obtain numerical values for the param- 
eters in the Hamiltonian. This makes it possible to ad- 
dress the band-to-acceptor transition for which there is 
no experimental data and no available prediction of how 
the PL spectra should look like. 



A. Band-to-band transitions: 4-spin-model 

In the band-to-band transition the system is left in 
one of the GSM states discussed in the previous section 
where the Mn-acceptor hole complex behaves like a F = 1 
spin. In band-to-band transitions the spin A complex is 
probed by the quantum dot exciton. This is the scenario 
considered by Kudclski et aL— . Here we compute the 
PL spectrum by numerical diagonalization of the GSM 
and the XSM within the 4-spin-model and then combin- 
ing equations ([3]) and the polarization operator equation 
(TT5|) . We use the numerical values of the coupling con- 
stants in the 4-spin-model as to fit the experimental PL 
spectrum of reference 0. The experimental zero field PL 
features 5 lines: a central low intensity one in between 
a high and a low energy doublet. Both the B z = zero 
field PL and the PL as a function of energy and applied 
field in the Faraday geometry are shown in figure [3J Our 
calculations, shown in figure ([3]), reproduce the zero field 
PL£ with 5 peaks at zero magnetic field, as well as the 
main features of the PL as a function of magnetic field. 
Notice that in the horizontal axis we plot with respect to 
Eq the transition energy of the bare quantum dot exciton, 
excluding its coupling to the Mn. 

The origin of the 5 peaks at zero field can be under- 
stood by inspection of the energy diagram shown in fig- 
ure (0) (sec also reference H), in which we only show 
states within the F = 1 manifold. The basic idea is 
that the quantum dot exciton is probing the F z compo- 
nent of a spin F = 1 object. If F z was a good quantum 
number, 3 lines should be expected: the middle peak, 
corresponding to F z = 0, and a high and low energy 
peaks, corresponding to the spin splitting of the F z = ±1 
states Ising coupled to the exciton spin. However, the 
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FIG. 3: Neutral exciton band-to-band transition calculated 
with the 4-spin-model. Upper panel: B z — a + -PL. Lower 
panel: color plot of a + -PL intensity as function of energy 
(vertical axis) and applied magnetic field (horizontal axis). 
(Color online) 

in-plane anisotropy results in the mixing states within 
the F z = ±1 doublet into x and y states with slightly 
different energies. As a result, there are direct (xx,yy) 
transitions as well as crossed transitions {xy,yx). The 
fact that the height of the xx and xy transitions are sim- 
ilar denotes that the exchange interaction is comparable 
to the in-plane anisotropy term E. The energy difference 
between these satellite (F z = ±1) transitions and the 
F z = central peak arises from the exchange coupling to 
the exciton. 

The origin of the reduced intensity of the F z = cen- 
tral peak is not in the quantum mechanical matrix ele- 
ments, but in the smaller statistical occupation proba- 
bility, given that the F z = state has a higher energy 
than the F z = ±1 doublet. Interestingly, since there are 
no transitions mixing F z = ±1 to F z = 0, D can not be 
inferred directly from PL line splittings. In contrast, the 
value of D stronlgy affects the intensity of the central 
peak. We estimate D ~ 4 meV. 

These results arc different from single Mn doped CdTe 
quantum dots, for which the PL has 6 peaks at zero field. 



There the QD exciton is probing the M z component of a 
spin 5/2 object without in-plane magnetic anisotropies. 
Here the 5 peaks show the interaction of a QD exci- 
ton with a spin F = 1 object with in-plane magnetic 
anisotropy. 

Additional information is obtained from the evolution 
of the PL spectra as a function of an applied magnetic 
field along the growth axis z. In the experiment^ the 5 
lines seen at zero field evolve, changing both in inten- 
sity and energy, in an intricatcd manner. In figure the 
lower panel of we show a contour plot of the PL inten- 
sity as a function of energy (vertical axis) and magnetic 
field (horizontal axis) for er + transitions, obtained within 
the 4-spin model. The fact that the calculation is in 
fairly good agreement with the experiment, and provides 
a strong back up for the theory. 



B. Band-to-band transitions: 2-spin-model 

We now discuss the physical interpretation of the evo- 
lution of the PL spectra as a function of the applied mag- 
netic field using the 2-spin-model proposed by Kudelsy et 
al&. The 2-spin-model affords analytical expressions for 
the PL spectrum at finite magnetic field in the Faraday 
configuration. The derivation is shown in the appendix. 
We address the merger of 3 lines at a particular value 
of the applied field, B*, the non-monotonic evolution of 
the highest and lowest energy lines at small field and the 
quenching of their intensity at large fields. Within this 
model a given circular polarization of the photon fixes 
the QD exciton spin. There are 3 exciton states and 3 
ground states. The F z = state, both for the XSM and 
the GSM, is decoupled from the other states and gives 
rise to the central line. The Zeeman shift of this line is 
that of the QD exciton, gx^BB z . Thus, we can fit the 
experimental data£ and infer from here the ^-factor of 
the QD exciton gx = 1.2 not far from values reported 
before^. 

The other four lines come from transitions within the 
F z = ±1 doublets. The energy levels of the ground state 
F z = ±1 doublet are given by —V ± he where 

h G = V^ 2 + (g F fi B B z ) 2 (22) 

The energy levels of the F z = ±1 doublet in the XSM 
with exciton spin X z = ±1 are given by Eq — T> + 
gx^BX z B z ± hx where Eq is the exciton energy transi- 
tion without spin and Zeeman terms and 

h x = ^£ 2 + {g F H B B z +JX z y (23) 

Both in the GSM and XSM the two states in the doublet 
are linear combination of both F z = +1 and F z = — 1. 
The mixing between the F z = ±1 states is governed by 
the competition between the in-plane anisotropy £, and 
the longitudinal interactions of the F z spin with the ap- 
plied field and, in the exciton manifold, the exchange 
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coupling to the exciton. This competition can be de- 
scribed by two angles, cot#c = 9FfJ -^ B = and cot9 x = 
9f^bBz+JX z Duc tQ thc diffcrcnt m i x i ng i n tne GSM 

and XSM, 4 transitions are allowed. We label them with 
ba, where a = ± and b = ± label the low (— ) and high 
(+) energy state of the ground and exciton state in the 
F z = ±1 manifold. From equation (|A13|) we have the 
transition energies 



E" 



Eo + 9x^bX z B z 



bh, 



ahr 



(24) 



At B z = we have h x = yj£ 2 + (JX Z ) 2 and h G = £, 
so that hx — he > 0. Thus, from the B z = point we 
immediately can label the 4 non-monotonic lines from low 

to high energy, at zero field, as — h (1), (2), ++ (3) 

and H — (4). We also denote with (0) the central line with 
F z = 0. The splitting between the two low (high) energy 
lines is 2ha and their average energy is — hx {+hx)- At 
zero field h G = £, thus £ is half the splitting within 
both the low and the high energy doublets. We thus 
infer £ = 0.035 meV. The splitting between the high 
energy doublet and the low energy doublet is, at zero 
field, 2h x = VJ 2 +£ 2 - From here we infer \ J\ = 0.14 
meV. 

Equation (f2"4"|) permits to extract the field B* at which 
lines (2), (3) and (0) cross. The crossing arises from the 
compensation between the Zeeman splitting of thc F spin 
and its exchange coupling to the exciton. The condition 
is hx = h G is satisfied for 



2g F /j, B B* 



-JX Z 



(25) 



Since B* is positive for X z = +1 we immediately see that 
J must be negative: the F z spin is fcrromagnctically cou- 
pled to the QD exciton. As discussed earlier, the negative 
sign can be obtained even if in thc 4-spin-model thc QD 
hole-Mn coupling is antiferromagnctic. Thus, the nega- 
tive sign comes from the strong correlation between the 
Mn spin and the acceptor hole spin, such that the sign of 
F z and M z are anticorrelated in the F = 1 manifold. In 
the simulations with the 4-spin model we have used pos- 
itive values for the QD holc-Mn coupling, obtaing good 
agreement with the experiment. Since we can infer J 
from the B z = data, equation (f25|) permits to infer gp 
from the experimental value of B* . We obtain gp = 3.01, 
not far from the g F = 2.77 of Mn-acceptor complex in 
GaAs. 

The intensity of the 4 lines in the F z = ±1 manifold 
are a function of a = h {9 G — 9 X ) (see equation IA13|) . 

In particular, the strength of H — \- and transitions 

is given by cos 2 (a) whereas the strength of thc H — and 
— h transitions is given by sin 2 (a). In the high field limit, 
when g F fiBB z 3> \ JX Z \, we have 6x = @G and a goes to 
zero. In figure (j4j we plot the a + PL spectra. The size 
and color of the circles is proportional to the quantum 
yield equation (|A13|) . The quantum yield of the F z = 
transition is constant. The slope of this line, coming from 
the Zeeman splitting of the QD exciton, is also present 
in thc other four lines. 
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FIG. 4: a -PL intensity of the neutral exciton band-to-band 
transition, as a function of the magnetic field, calculated with 
the 2-spin-model (eq. IA12|) . Thermal effects are not included. 
The size and color of the symbols of the numbered lines is pro- 
portional to the quantum yield. The intensity of the central 
line is constant. (Color online) 



Lines (1) and (4) come from transitions that mix + 
and — states with different symmetry. Their energy with 
respect to line (0) is given by ±(hx + ha). Since J < 0, 
hx has a minumum at B z = B* whereas he, whose con- 
tribution is smaller, has a minimum at B z = 0. The 
intensity of these lines quenches as the magnetic field in- 
creases so much that g F [iBB z S> £ and the F z is restored 
as a good quantum number. In contrast, lines (2) and 

(3) come from ++ and transitions. At large fields 

the quantum yield is increased since the mixing between 
F z = +1 and —1 is quenched. The energy of (2) and (3) 
with respect to (0) is given by ±(hx — ho). Thus, when 
the Zeeman splitting is much larger than the exchange 
coupling, i.e. for B z ^> 2B* , we have hx — h G and the 
slope of lines (2) and (3) is the same as (0). 

The model captures the main experimental features of 
the PL spectrum^, namely: i) 5 peaks distributed as a 
high and low energy doublets with a small intensity cen- 
tral peak; ii) as a magnetic is applied along the growth di- 
rection, the central line does not change intensity and has 
a linear shift, whereas the doublets have non-monotonic 
shifts and do change intensity, two of them fade away; 
Hi) at a given value of B* 3 lines, coming from the low 
and high doublets and the central line, cross. 

The results of the 2-spin-modcl (figure |4|) and 4-spin- 
model (figure [3]) have no apparent differences (besides 
the lack of thermal occupation in the 2-spin-model). This 
validates thc approximations made to go from the 4-spin- 
model to the 2-spin-modcl. 

The model portraits the neutral Mn acceptor complex 
in InAs as a spin F = 1 nanomagnet with two almost 
degenarate ground states, F z = ±1, a rather large sin- 
gle ion magnetic anisotropy of X>, and a small in-plane 
anisotropy £ . The quantum dot exciton is Ising coupled 
to F z and permits a direct measurement of £ and 
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and indirect measurement of T>. Finally, we have verified 
that the model proposed by Climente et al£^ could not 
account for the experimental PL reported by Kudclski 
et ai. In a nutshell, this model is very close to the one 
proposed by one of us to account for the PL of Mn doped 
charged CdTe quantum dots^. In the model of Climente 
et ai, the ground state manifold has 6 doublets com- 
ing from the Ising coupling of the hole and the Mn spin. 
The injection of an additional electron-hole pair will re- 
sult in two holes with opposite spin occupying the same 
orbital, uncoupled from the Mn, which would interact 
only with the electron, presumably via a Heisenberg cou- 
pling. Thus, the exciton manifold of such a model would 
have two spectral lines. The resulting spectrum would 
have 11 lines with a characteristic V shaped. If the the 
Mn-clcctron coupling is turned off, the model yields 6 
equally strong lines. 



C. Band-to-acceptor transition 

We now consider band-to-acceptor transitions, 
X°A° — > h + A~ such that, after photon emission there 
is a hole left in the QD levels and the Mn is liberated 
from the acceptor hole. This kind of transition has 
been observed in Mn doped GaAs bulb^ and quantum 
wells^ 3 . but not yet in Mn doped quantum dots. Three 
obvious differences with the band-to-band transitions 
can be mentioned beforehand. Firstly, the PL spec- 
trum associated to the acceptor transition should be 
red-shifted with respect to the band-to-band transition, 
by the sum of the acceptor binding energy and the 
quantum dot confinement energy. In GaMnAs quantum 
wells the reported shift is approximately 107 meV— . 
Secondly, we expect a smaller intrinsic efficiency of the 
band-to-acceptor process compared to the band-to-band 
case, due to the smaller electron-hole overlap in the 
case of the former. Thirdly, the final state, A~h + is an 
excited state since the QD hole could be promoted to 
the acceptor state, reducing the energy of the system. 
The spin of the A~h + state is the product of the Mn 
spin 5 = 5/2 and the QD hole spin. 

A diagram of the energy levels of this transition is 
shown in figure The states of the XSM are the same 
as in the band-to-band transition, except for the fact that 
now both the X = ±1 and the X = ±2 transitions are 
allowed, since the optical selection rule must be enforced 
with the acceptor hole and not with the QD hole. Quan- 
tum dot electron-hole exchange splits the "dark" and 
"bright" lines. Thus, in the lowest energy F = 1 mani- 
fold there are 6 energy lines in the XSM corresponding 
to 3 projections of F z and the two excitons. 

The GSM features now an ionized Mn acceptor inter- 
acting with a quantum dot hole instead of with an ac- 
ceptor hole. Both the strength and the symmetry of this 
coupling arc different: the QD hole-Mn coupling is much 
weaker than acceptor hole-Mn coupling and, due to the 
lack of spherical symmetry of the QD hole state, is pre- 




wcak lines due to mixing by £ [F£ — F?) 



FIG. 5: Energy level diagram for the neutral exciton band-to- 
acceptor transitions. Both the ground state and exiton state 
energy lines are shown, as well as their evolution as a function 
of the magnetic field. (Color online) 

dominantly Ising^^. Hence, the GSM is given by the 
Ising coupling of the ionized Mn, with spin S = 5/2 and 
the QD hole, with total angular momentum = ±3/2. 
The spectrum of this system, relevant for single Mn in 
CdTe QD 11 ' 13 , is formed by six doublets. We can label 
the GSM states with the projections of the Mn spin and 
the QD hole spin along the growth axis, \M z ,ah)- Their 
eigenvalues are: 

E Mz .a h =+J^ M M z a h . (26) 

where J^ M is the exchange coupling between the QD hole 
and the ionized Mn spin. 

There could be as many as 36 PL spectral lines joining 
the 6 energy levels of the XSM and the 6 energy levels 
in the GSM. The highest energy PL would correspond 
to the highest energy exciton state (within the F = 1 
manifold), with quantum numbers F z = 0, X = ±1 and 
the lowest energy ground state, with quantum numbers 
M z = +5/2, a h = -3/2 or M z = -5/2,a h = +3/2. 
Interestingly, the M z = ±5/2 state has zero overlap with 
the F z = 0, so that this particular transition is forbidden 
(see figure[S]). Going up in the M z ladder, the first excited 
states in the GSM are M z = +3/2, a h = -1/2 and M z = 
-3/2,(7,, = +1/2. 

The results of the simulation within the 4-spin-model 
are shown in figure ([6]). Using values from the band to 
band transition we take D = 4.7 meV, E = 0.31 meV, 
e = 5 meV, J hM = -0.0405 meV, J eh = -0.2 meV. 
We assume that the QD hole-acceptor hole coupling is 
zero, but we take a ferromagnetic coupling between the 
Mn and the QD hole to reproduce the effective ferromag- 
netic coupling between the exciton and the Mn acceptor 
complex. In the ground state we take a QD hole-Mn 
coupling larger than that of the XSM to account for the 
larger electrostatic attraction of the ionized Mn acceptor, 
JhM = 0-15 meV. The PL features a group of higher inte- 
sity lines at low energy and a group of weaker lines 2 meV 
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above. As we discuss below, the splitting between the 
low and the high energy groups turns out to be given by 
T> = 3D/10. Thus, band-to-acceptor transitions would 
permit a direct spectroscopic measurement of D. The V 
shape of the intensity pattern in the low energy group is 
somewhat related to the PL spectra of charged Mn-doped 
CdTe quantum dotst In both cases the optical matrix 
elements feature overlap between states with well defined 
Mn spin and states where the Mn spin is Heisenberg ex- 
changed coupled to another carrier (the extra electron in 
Mn doped CdTe and the acceptor hole in the case con- 
sidered here). 
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FIG. 6: B z = a -PL for neutral exciton band-to-acceptor 
transition, as calculated with the 4-spin-model. (Color online) 



D. Band-to-acceptor transition: 2-spin-model 

The results of figure (O can be rationalized with the 
2-spin-model, like in the band-to-band transition. We 
describe exciton states as the product of the QD electron- 
hole pair spin, which can be ±2 or ±1 times the acceptor 
complex spin F z = ±1,0. We ignore, for a moment, the 
in-plane anisotropy term £ that mixes F z = ±1. Thus 
F z , a c and ah are good quantum numbers in the XSM 
and M z and ah are good quantum numbers in the GSM. 
Of course, the spin of the QD hole, ah, and the Mn spin, 
M , are conserved during the transition. 

The quantum matrix elements of these transitions are 
given by the matrix elements of the dipole operator (eq. 
[20land[2~lj) between emitting states, X and ground states 
with quantum numbers M Zl a%: 

\(M z ,a h \V + \^ x )\ 2 = \C^ t+3/2tl , ah + ^C^, +1/2 , T ,J 2 

\(M z ,a h \r^ x )\ 2 = |cf /2i _ 3AT , CTh + -^,_ 1/2 , i)CT J 2 

for the a + and a~ transitions respectively. Ignoring the 
in-plane mixing term £ , for a given QD exciton (a c , ah), 



the coefficients ^ a a arc given by equation (|5|). 
Thus, as long as F z is a good quantum number, the spin 
of the acceptor hole and the Mn satisfies the rule 

M z =F Z - Jz (27) 

Thus, for a given state in the XSM, with quantum 
numbers \F = 1, F z , a c , ah), and a given polarization of 
the photon, ±1, we immediately get the permitted values 
of the annihilated acceptor hole spin j z = ±1 — a c (eq. 
[HI)) and the allowed values of the Mn spin after photon 
emission (cq. |2T|) . 

Using equations (J5]), (fl"9|) . ([27)) one can make a ta- 
ble where the inputs are the QD electron spin a c , the 
acceptor complex spin F z , the circular polarization of 
the photon a^, and the outputs are the annihilated ac- 
ceptor hole spin j z , the Mn spin M z after photon re- 
combination and the quantum yield of the transition, 
X = I(F Z , a c ,a^, M z ). For a + polarization we obtain: 



F z 


a c 


jz 


M z 


T(<7+) 


+1 


T 


+ 1/2 


+ 1/2 


5/100 





T 


+ 1/2 


-1/2 


10/100 


-1 


T 


+ 1/2 


-3/2 


10/100 


+1 


1 


+3/2 


-1/2 


5/100 





1 


+3/2 


-3/2 


20/100 


-1 


1 


+3/2 


-5/2 


50/100 



The table for a~ is obtained by application of the 
time reversal operator, which changes F z ,a c ,j z ,M z , to 
—F z ,—a c , —j z , —M z , and give the same intensity. In or- 
der to get the PL spectrum we need the transition ener- 
gies Ex — Eg- Neglecting the £ term that mixes F z = ±1, 
the XSM energies are approximated by 

Ex (F z ,a c ,a h ) = E - VF 2 Z + JX Z F Z + J eh a c a h (29) 

and the energy of the ground states by equation (|26p . 
Thus, the transition energies uj = Ex — Eq are given by: 

u = Ev-VF 2 +J{a c +ah)F z +ah{Jehac-Jh M M z ) (30) 

In figure ([7|) we plot the intensity of the lines, with- 
out thermal depletion of the higher energy F z = states. 
The numerical values of the constants are those of the 
4-spin-model simulation, properly rcnormalized to the 2- 
spin-model. Thus we take V = 3D/10 = 1.41 meV, 
j = 21.7,, aj = _Q.12 me v. The numerical values of the 

o 

QD electron-hole exchange and the QD hole-Mn coupling 
in the ground state are the same in the 2-spin and 4-spin 
models, J eh = -0.2 meV and J,f M = +0.15mcV. Wc 
take Jf^ M > J because the quantum dot hole is electro- 
statically attracted towards the Mn when this is ionized. 
The PL of figure © feat ures 12 peaks, corresponding to 
the 6 allowed transitions of equation (f2"5| times the two 
possible spin orientations of the QD hole. The high en- 
ergy group corresponds to the transitions coming from 
the F z = states. Most of the splitting between the high 
and the low energy group is given by T>, with some contri- 
bution coming from the exciton Mn-complex exchange. 
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The lower energy group of 8 peaks correspond to transi- 
tions in the F z = ±1 manifold. The high intensity peaks 
correspond to transitions where the Mn spin is left in a 
±5/2 state. For <j + transitions, shown in the figure, the 
peaks correspond to final states with M z = — 5/2 and the 
quantum dot hole with spin +1/2 (low energy peak) and 
— 1/2 (high energy peak). The splitting between those 
two is 5 l M + ~ 3 • Interestingly, this simple model 
captures the main features of the PL, as calculated with 
the 4-spin-model. The latter has more lines because of 
the the 8 terms, that mix F z = 1 to the F z = — 1 terms. 
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FIG. 7: B z = <j + -PL for neutral exciton band-to-acceptor 
transition, calculated with the 2-spin-model with £ = (see 
text). (Color online) 



It is worth noting that the band-to-acceptor transi- 
tions permit to relate the degree of spontaneous cir- 
cular polarization to a possible inbalancc in the F z 
populatio n 33 i 34 i 35 . For instance, the transitions that have 
has initial state F z = —1, averaged over all the possible 
QD spin orientation and final state have a degree of cic- 
ular polarization of more than 70 percent: 



£ Mg , c W-l. g o M z ,+) - /(-!, <7 C , M z , -)) = 5 
Em, a c (H-l,<rc, M z ,+) + a c , M z , -)) ~ 7 

(31) 

Thus the degree of Mn spin polarization can be probed 
by measuring the degree of circular polarization of the 
PL. 

Hence, whereas band-to-band transitions probe the 
Mn-acceptor complex, which behaves as a spin F = 1 
object, the band-to-acceptor transitions connect initial 
states for which the Mn spin is correlated to the accep- 
tor hole to final states where the Mn spin is Ising coupled 
to the QD hole, but with M z as a good quantum num- 
ber. In band-to-acceptor transitions the photon energy 
and polarization carry information about M z . 



IV. CHARGED EXCITON SPECTROSCOPY 



We now consider the PL spectra of negatively charged 
InAs dots. From the theory side there is a lot of inter- 
est on the effect of number of carriers on the magnetic 
properties of a dot doped with Mn atoms^ i 11 ' 15 ^ 4 ' 36 ^ 7 ' 38 . 
A priori, the ground state of a negatively charged InAs 
dot doped with a single Mn should be the ionized A~ 
acceptor, with spin properties identical to those of Mn in 
neutral CdTe 1 ^. Band-to-band transitions should yield 
zero field PL spectra with 6 peaks. Band-to- acceptor 
transitions for the negatively charged dot, e~A° — > A~ , 
have been studied by by A. Govorov in referenced^. He 
obtained a PL spectrum with 3 peaks in at zero magnetic 
field. 




A- 



FIG. 8: Possible electronic configurations for two photon de- 
cay of the negatively charged InAs QD with one Mn. The up- 
permost diagram shows the configuration with two QD elec- 
trons, one QD hole and the neutral Mn-acceptor hole complex. 
By a band-to-acceptor transition (pathway 1), one arrives at 
one of the intermediate configurations, with an ionized Mn 
and a QD exciton. A further band-to-band transition, cre- 
ates the lowest energy configuration: an ionized Mn. Another 
recombination possibility is via pathway 2: a band-to-band 
transition yields the intermediate configuration with the neu- 
tral Mn-acceptor complex and one QD electron. A further 
band-to-acceptor transition then yields the same end config- 
uration as the one achieved via pathway 1. (Color online) 
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Contrary to these expectations, the experimental re- 
sults on negatively charged single Mn doped InAs QD^ 
show very similar neutral and charged exciton band-to- 
band spectra, with 5 peaks in the zero field PL. The re- 
ported transitions correspond to an emitting state with a 
neutral Mn acceptor plus a quantum dot trion (X~A a ). 
Thus, there are two electrons in the conduction level, and 
two holes, the acceptor hole and the QD hole (sec figure 
[8]). The final state, after the emission of two photons, is 
the ionized Mn, with spin S = 5/2. The fact that the 
reported neutral and negatively charged transitions are 
very similar indicates that QD electrons are very weakly 
coupled to the Mn-acceptor hole complex. This is dif- 
ferent from CdTe dots doped with one Mn, where the 
addition of a single electron changes the spin properties 
of the Mnl ' 11 ' 15 . In InAs, the Mn is strongly coupled to 
the acceptor hole and is less sensitive to the number of 
conduction band carriers. 

Yet the presence of additional QD electrons can result 
in new PL spectra when we consider band-to-acccptor 
transitions, unreported so far in single Mn-doped InAs 
quantum dots. Since the negatively charged trion X~ A 
features 2 electrons and 2 holes in different states, the 
decay towards the ground state A~ can occur via two 
recombination pathways, as show in figure [5J One of 
the pathways (denoted as 1 in the figure) consist in a 
band-to-acceptor transition X~A° — > X°A~ followed by 
a X°A" — > A~ , i.e., a band-to-band transition that an- 
nihilates a QD exciton coupled to a ionized Mn acceptor. 
The PL of the first step in pathway 1 is related to the 
band-to-acceptor transition of the neutral dot discussed 
in the previous section. The PL of the second step in 
pathway 1 should be very similar to that of single Mn 
doped CdTe QD. As long as the Ising part of the cou- 
pling between the QD hole and the Mn spin is dominant, 
as in CdTe, the PL of the energy and polarization of the 
photon yield direct information of the Mn spin after the 
exciton recombination^. 

Pathway 2 in the figure starts with a band-to-band 
transition X~A° — ► e~ A followed by a band-to-acceptor 
transition e~A° — ► A~. The first step has been observed 
experimentally and, since the electron-Mn exchange is 
very weak, yields a PL very similar to the neutral case 
considered above. The second step is identical to the 
transition considered by Govorov^. 

In the upper panel of figure ^ we show the zero 
field PL for the band-to-acceptor transition for negatively 
charged InAs QD doped with 1 Mn. The emitting state is 
X~A° and the final state is X°A~ , a quantum dot exci- 
ton coupled to an ionized Mn acceptor. The correspond- 
ing energy level diagram is very similar to that of the 
neutral band-to acceptor transition shown in figure |5|). 
Ignoring states with F > 1, there are 6 exciton states, 
corresponding to 3 F z values and 2 QD hole spin states. 
There are 24 ground states, corresponding to the 6 Mn 
spin orientations and the 4 spin states of the quantum 
dot exciton. This is in contrast with the X°A° — > h + A~ 
transition, for which both the GSM and the XSM have 12 



0.004 



o 







" 0.050 








- 






0.035 " 










, A 


j i L . 






1. 





-35 -34.5 -34 -33.5 -33 -32.5 -32 -31.5 
E-E [meV] 




Magnetic field in z-direction [T] 

FIG. 9: Upper panel: B z = cr + -PL for a negatively charged 
exciton band-to-acceptor transition (X~ A — » X°A~), as cal- 
culated with the 4-spin-model. Lower panel: Intensity plot 
for the <t + -PL for the negatively charged exciton band-to- 
acceptor transition (X~A° — > X°A~), as a function of energy 
(vertical axis) and applied magnetic field (horizontal axis). 
(Color online) 



states. Other differences with that transition is the lack 
of QD electron-hole exchange in the emitting state, and 
the presence of that coupling in the ground state. Thus, 
electron-hole exchange is present both in the neutral and 
in the charged case, either in the XSM or in the GSM. 
It would be possible to do an analytical model for the 
charged case along the lines of the previous section for 
which the optical matrix elements would be still given by 
equation (gHJ). Not surprisingly, the PL X~A° — > X°A~ 
of figure ([9]), calculated with the 4-spin-model, is quite 
similar to that of the X°A° -> h + A~. 

In the lower panel of figure ([9]) we plot the er + PL as 
a function of the magnetic field in the Faraday config- 
uration. The two low energy brighter peaks correspond 
to transitions where the Mn spin is, in the final state 
M z = —5/2. They are splitted due to the different spin 
orientations of the QD exciton to which they are coupled. 
This results also in different slopes as the magnetic field is 
ramped. The evolution of the energy levels in the ground 



13 



state results in a compensation of the zero field exchange 
splittings by the finite field Zeeman splittings. At a par- 
ticular value of the magnetic field, several lines become 
degenerate. In the absence of spin-flip terms, they do not 
anti-cross. This particular energy arrangement has been 
reported in CdTe doped with Mn&. The anti-crossings 
observed at B z ~ IT are related to those taking place 
at the XSM and dicussed above in the context of neutral 
band-to-band transitions. 



V. DISCUSSION AND CONCLUSIONS 

We have addressed the problem of single exciton spec- 
troscopy of a single-Mn InAs doped quantum dot. The 
main goal is to link a few-spin Hamiltonian with the PL 
spectra featuring spin-split peaks at zero magnetic field. 
We have focused on the fact that, for single Mn-doped 
InAs QD this is a 4 body problem, with the QD elec- 
tron, QD hole, acceptor hole and Mn spin. In order to 
account for the experimental observations^ it is impor- 
tant to assume that a bulk-like acceptor states surives 
inside the gap, weakly affected by the quantum dot. The 
strongest exchange interaction is that of the Mn and the 
acceptor hole. In most instances this permits to interpret 
the results as if the quantum dot exciton interacts with a 
spin F = 1 object, obtained from the antiferromagnetic 
coupling of the Mn spin S = 5/2 and the acceptor hole 
j = 3/2. We use both a 4-spin-model, in which the iden- 
tity of the spin of all the carriers is included in the cal- 
culation, and a 2-spin-model, that ignores the composite 
nature of the exciton and the Mn-acccptor complex. The 
models are diagonalized numerically and the PL spectra 
are obtained, taking full account of the optical and spin 
selection rules. 

The 2-spin-model portraits the Mn acceptor complex 
in InAs as a spin F = 1 nanomagnet with two almost 
degenarate ground states, F z = ±1, a rather large sin- 
gle ion magnetic anisotropy of 2?, and a small in-plane 
anisotropy £ . The notion that single Mn atoms in quan- 
tum dots can behave like aritifical single molecule mag- 
nets has been discussed befor o 11 ' 39 . Interestingly, the 
zero field exciton spectroscopy gives a direct measure- 
ment of £ ~ 0.035 meV, but only an indirect measure- 
ment of D ~ 4 meV, through the height of the central 
peak. The evolution of the PL spectra as a magnetic 
field is applied permits an indirect measurement of other 
energy scales in the problem. In the band-to-band tran- 
sitions, measured by Kudelski et al4, there is particular 
value of the field B* at which three lines in the spectra 
merge. B* is the field at which the Zeeman splitting and 
exchange coupling to the exciton have the same intensity 
and opposite sign (eq. I25[) . Since the zero field measure- 
ment provides 3 , the value of g F in InAs can be inferred 
from B*. We estimate g F = 3.01. 

The nature of Mn spin S — 5/2 can be unveiled in two 
manners: upon electron doping the system and in band- 
to-acceptor recombination of the charge neutral dot. The 



latter results in the ionization of the Mn complex so that, 
in the final state, the Mn spin is S = 5/2 and is presum- 
ably Ising coupled to the QD hole. The band-to-band 
and band-to-acceptor transitions are very different. In 
the former the Mn spin is slaved by the acceptor hole both 
in the XSM and GSM states. Thus, photon emission does 
not change the symmetry of the Mn spin Hamiltonian. 
In the band-to-acceptor transition the Mn is liberated 
from the acceptor hole in the final state, so that pho- 
ton emission involves a change of the effective Mn spin 
Hamiltonian. In this sense, the band-to-acceptor tran- 
sition resembles the negatively charged trion transition 
in Mn doped CdTei. Band-to-acceptor transitions would 
provide a direct spectroscopic measurement of T>. 

Notice that, within the 2-spin-model we take an effec- 
tive ferromagnetic coupling between the exciton and the 
Mn. The 2-spin-model does not say whether this coupling 
is QD hole-Mn acceptor complex, QD hole-acceptor hole, 
QD clcctron-Mn or QD electron-acceptor hole. If we as- 
sume that the dominant coupling is between the Mn spin 
and the QD hole, we could conclude that the coupling is 
ferromagnetic, at odds with the usual antiferromagnetic 
coupling of holes and Mn in III-V materials. Interest- 
ingly, we have seen how the bare sign of the QD hole-Mn 
coupling is reversed when going from the 4-spin-modcl 
to the 2-spin-model (see equation [TT)) . The strong Mn- 
acceptor hole complex results in the renormalization of 
the spin interactions of these two spins with applied field 
and the exciton spins. This can even result in sign inver- 
sion of the exchange interaction. Thus, we have shown 
that a ferromagnetic coupling between the 2 spins of the 
simpler model can arise even if the underlying spin cou- 
plings of the QD hole to the Mn spin are antiferromag- 
netic. In this way, we reconcile the standard view about 
this system, in which the holes are coupled antifcrromag- 
netically to the Mn spin, with the the observations^. 

A single Mn in a quantum dot provides an ideal system 
to address and control the spin of a single object in a solid 
state environement. The PL spectra provide valuable in- 
formation of the effective spin Hamiltonians, for different 
states of the dot. A Mn atom in a InAs QD behaves like 
a 3 level system. It might be possible to use the encode 
a single qubit in the almost degenerate F z = ±1 doublet. 
The presence of higher energy F z = state, acting as a 
barrier, might block spin relaxation between the ground 
states. In analogy with the case of electron doped and 
hole doped quantum dots, it should be possible to manip- 
ulate the spin of a single or a few Mn spins in a quantum 
dot by application of laser pulses. 
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APPENDIX A: ANALYTICAL SOLUTION OF 
THE 2-SPIN-MODEL 

Here we provide an analytical solution of the ground 
state and exciton state manifolds within the 2-spin- 
modcl. The GSM has 3 states, but the Hamiltonian can 
be block diagonalizcd since the F z — state is decou- 
pled from the F z = ±1 doublet. The XSM state has 
12 states, corresponding to the 4 spin orientations of the 
exciton and the 3 states of F z in the F = 1 manifold. 
Since we consider Ising interaction between the exciton 
and the F spin, the Hamiltonian of the XSM is also block 
diagonalizcd. Importantly, both the GSM and the XSM 
are each characterized by a single angle, G and Ox i that 
characterize the ratio between the in-plane mixing of the 
F z = ±1 components and their splitting, induced both 
by exchange interaction with the exciton and by Zeeman 
coupling. As we show here, the lineshape of the PL spec- 
trum depends on | (fix — Og)- 

In the basis |1, +1), |1, — 1), |1,0) the Hamiltonian of 
the GSM reads: 



H G = 



-V + c]fHbB z 
£ 




£ 
-V-g Ff x B B z ) (Al) 




We can write the two by two matrix within the F z = ±1 
subspacc as a linear combination of the unit matrix J, 
and the Pauli matrices <x z and er x : 

H G> ±i = -VI + h G (cos(6 g )ct z + sm(0 G )(T x ) (A2) 

with 

h G = (£, 9f^bB z ) = h G (sin(0 G ), cos{0 G j) (A3) 
where h G = £ 2 + (gfj,BB z ) 2 . Thus, we have 



,/ <, 9f^bB z 
cot(0 G ) = - 



(A4) 



At zero field we have 6* G = ir/2. 

The eigenvalues of the GSM are: —V — h G , —T> + h G 
and 0. The corresponding eigenvectors are the product 
the quantum dot ground state, denoted by |0), and the 
spin part, denoted by ^jfc and ^ G respectively. The 
spin part of the F z = ±1 sector reads 



* G 



sin(f; 



cos 



-cos(^) 
sin(%) 



(A5) 



Notice that, in the limit of very strong field, g F [iBB z 
£ , the mixing between F z = +1 and —1 vanishes. 

We now consider the XSM, which is split in 4 sec- 
tors with 3 states each and a well defined exciton state 
X z = ±1 and X z = ±2. We focus on the optically ac- 
tive excitons, X z = ±1. Since X z commutes with the 
XSM Hamiltonian, the X z = — 1 and X z = +1 sectors 
decouple and are described by 3 by 3 matrices with the 
same structure as (|A1|) . The F z — <E> X z states are 



decoupled from the F z = ±1 <g> X z states. In the basis 
|1,+1) ® |A 2 ), |1,-1) ® \X Z ) the XSM Hamiltonian for 
exciton X z reads: 

H Xz =±i = (E{X Z )-V)I 

+h x (cos(6 x )<T z + sm(6 x )cT x ) (A6) 

where Ex(X z ) = E a + gx^BX z B z is the energy of the 
excitonic transition neglecting spin coulings plus the QD 
exciton Zeeman splitting, and 



h x = (£,g F ^BB z + JX Z ) 
= h x (sm(6>x),cos(6»x)) 



where 



h x = yj£ 2 + (g F ^BB z + JX z y 



Now we have 



cot(0 x ) 



9fHbB z + JX Z 
£ 



(A7) 



(A8) 



(A9) 



Notice that the angle Ox depends both on the magnetic 
field B z and on the exciton spin projection, X z = ±1. 
At zero field the angles corresponding to X z = +1 and 
X z = — 1 differ by tt, so that cot($x) = At finite 

fields the effect of exchange and magnetic field on the Mn- 
acceptor hole complex can either compete or cooperate 
and the relation between 9x for X z = ±1 is non-trivial. 
The eigenvalues are for the XSM are: 



Eq 

Eo 



gx^BX z B z 
gxHBX z B z 
9x^bX z B z 



V 
V 



h x 
hx 



(A10) 



The corresponding eigenvectors are the product of the 
quantum dot exciton, denoted by \Xz) and the spin part, 
denoted by ^x> ^x respectively The spin part of 
the XSM eigenvectors in the F z = ±1 sector are: 



x 



sin ( 

cos I 



-cos(^) 
sin " 



X J \ VJ "° \ 2 J V 2 

The intensity of the PL lines is given by 



E 

a . b 



n b \(n\^x)\ 2 S(E x 



E a G 



(All) 



(A12) 



where both a and b run over +1, —1, 0. Here rib is the sta- 
tistical occupation of the exciton states. Both the matrix 
elements and the allowed transitions depend on the exci- 
ton spin X z = ±1, which is in turn given by the photon 
polarization. The intensity table | {^ G \^x) | 2 reads: 







+ 







cos 2 




G - Ox)) 


sin 2 


:k 


G - Ox)) 





+ 


sin 2 


'M 


G - Ox)) 


cos 2 




G - Ox)) 
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This matrix has 5 non-zero elements, corresponding to 
the 5 permited transitions for a given circular polariza- 
tion. The transition energies within the F z = ±1 dou- 
blets are given by 

E b ^ a = E + g x ^ B X z B z + bh x - ah G (A13) 

where both a and b can take values equal to ±1. 

As expected, the state is decoupled from the others 
and the optical matrix element is 1. The 4 transitions on 



the F z = ±1 doublet ( , — h, + — ,++) are governed 

by a single variable \ (9q — Ox)- The intensity of the 
crossed sign transitions goes to zero in the limit of very 
large fields or exchange coupling. In this situation the 
mixing induced by in-plane anisotropy term, E, is neg- 
ligible. Notice that the height of the actual transitions 
depends both on K^gI^^)! 2 and on the statistical occu- 
pation. Thus, the observed intcntisty of the F z — line 
is smaller, due to a smaller statistical occupation. 
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